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We study theoretically the spin dynamics of antiferromagnetic molecular rings, such as the ferric 
wheel Feio- For a single nuclear or impurity spin coupled to one of the electron spins of the ring, we 
calculate nuclear and electronic spin correlation functions and show that nuclear magnetic resonance 
(NMR) and electron spin resonance (ESR) techniques can be used to detect coherent tunneling of the 
Neel vector in these rings. The location of the NMR/ESR resonances gives the tunnel splitting and 
its linewidth an upper bound on the decoherence rate of the electron spin dynamics. We illustrate 
the experimental feasibility of our proposal with estimates for Feio molecules. 
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Magnetic molecular clusters have been the subject of 
intense research in recent years, because they offer the 
possibility of observing macroscopic quantum phenom- 
ena |fflj. Ring systems such as the antiferromagnetic 
(AF) ferric wheels, Fe6 and Feio l@~fi' a ^ so allow one 
to study the transition from microscopic magnetism to 
one-dimensional bulk magnetism. In ferromagnetic clus- 
ters such as Fes and Mni2, incoherent tunneling of the 
S = 10 magnetic moment is observed directly in mag- 
netization and susceptibility measurements (|||. The 
AF ferric wheels are candidates for the observation of 
macroscopic quantum coherence (MQC) in the form of 
coherent tunneling of the Neel vector |Q . Although quan- 
tum effects in antiferromagnets are expected to be more 
pronounced than in ferromagnets @||, the detection of 
quantum behavior is experimentally more challenging . 
The reason for this is that magnetization and suscep- 
tibility measurements probe only the total spin of the 
molecule which, by symmetry, remains unaffected upon 
tunneling of the Neel vector n. At low temperatures, the 
dynamics of n is determined by two characteristic fre- 
quencies, the tunneling rate A/h and the electron spin 
decoherence rate Ts- In Feio, A can be tuned from 
to 2K by varying the magnetic field @]. Although the 
tunnel splitting A enters thermodynamic quantities such 
as magnetization and specific heat, no conclusive results 
revealing the characteristic functional dependence A(£?) 
have yet been obtained || . Measuring Ts is of central im- 
portance for the characterization of the quantum tunnel- 
ing, since coherent quantum tunneling requires < A. 
An estimate of Ts can be obtained from the typical en- 
ergy scales of various interactions leading to decoherence. 
Spin-phonon interactions are frozen out at low T. Nu- 
clear dipolar (O.lmK) and hyperfme (lmK) interactions 
are significantly smaller than interring electron spin dipo- 
lar interactions (some lOmK). These numbers indicate 
that, in Feio, tunneling of n is coherent for a wide range 
of A. 

Motivated by these numbers, in this paper we study 



the quantum spin dynamics of the ferric wheel. We show 
that the tunnel dynamics of n enters the correlation func- 
tions of a single electron spin. As indicated in |i"o[ for a 
bulk AF system, a nuclear spin coupled to a single elec- 
tron spin acts as a local spin probe. We show that in 
small AF rings in which n itself has additional coherent 
dynamics, the correlation functions of the nuclear spin 
exhibit signatures of the tunneling of n. In particular, 
we discuss the coherent dynamics of one nuclear spin cou- 
pled to one of the electron spins of the ferric wheel and 
show that both the tunnel splitting A and the electron 
spin decoherence rate can be obtained from NMR and 
ESR spectra. 

The ferric wheels Fe6 and Feio are well character- 
ized glllnjQ. The s = 5/2 Fe III ions are arranged 
on a ring, with an AF nearest-neighbor exchange cou- 
pling J > and a weak, easy-axis anisotropy directed 
along the ring axis e z . The minimal Hamiltonian for the 
system is 
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where N — 10 or 6 and sjv + i = si, h = g/isB, with B 
the external magnetic field and 3 = 2 the electron spin 
(/-factor. For Feio, J = 15.56K and k z = 0.0088J. For 
Fe6, the values for J and k z vary appreciably depending 
on the central alkali metal atom: for Na:Fe@, J = 32.77K 
and k z = 0.0136J, whereas for Li:Feg, J — 20.83K and 
k z = 0.0053J 

For k z — 0, the eigenstates of Hq are also eigen- 
states of the operator of total spin, S = X)iLi s ii with 
E s ,s m = (2J/N)S(S + 1) + hS x §. For k z < 2J/(iVs) 2 , 
the anisotropy can be taken into account in perturba- 
tion theory in k z . The scenario changes drastically for 
k z > 2J/(Ns) 2 , when n (staggered magnetization) is 
localized along ±e z . We label the states with n ori- 
ented along ±e z by | f) and | J.). In a semiclassical 
description, the low-energy sector of the ferric wheel con- 
sists of two states, a ground state, \g) = (| t) + I J- 
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))/V2, and a first excited state, |e) = (| T) - | 4))/V2. 
The static equilibrium properties of a system described 
by Eq. (Q) are discussed in detail in Ref. With 
B || e x , i.e. in the plane of the ring, the system 
exhibits interesting spin dynamics: In the high field 
regime, h x 3> Hloq, n is confined to the (y, z)-plane and 
tunneling takes place through the potential barrier of 
height Nk z s 2 . In particular, A = Ao| sin(7riV/i ;c /4J)|, 
with Ao = 8fiwo e~ s / h , lhq = s^/8Jk z /h, and 
S/h = Ns\/2k z / J, shows oscillatory behavior as a func- 
tion of h x 0- We restrict our considerations henceforth 
to the geometry B || e x . For Feio, Ao — 2.18K is much 
larger than in, e.g., Mni2 HJJlJI or Fes- 

Spin susceptibilities. - We consider first the stan- 
dard ac. spin susceptibility and ESR measurements, 
in which an infinitesimal magnetic probing field cou- 
ples to the total spin S = J2i=i s » 01 the ferric wheel. 
We will show that these experimental techniques are 
insufficient to detect coherent tunneling of n in a sys- 
tem described by Hq alone. In an effective-action de- 
scription for the system (|lj) with Lagrangian density 
L E [n] 0, we find S = -j^[mxn-h + n(n • h)]. In 
high magnetic fields h x 3> fcuiQ, the spin susceptibil- 
ity, Xaa{u) = Jo dr (T T S a (T)S a )e tulnT \i Un ^ u+i o, with 
a = x,y,z, (3 = l/fc^T, and T T the imaginary time order- 
ing operator, may be evaluated using spin path integrals. 
For uj,kBT/h<^ ujQ,h x /h, up to corrections 0(e~ s / h ), 
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probes allow the observation of the Neel vector dynamics. 
Nuclear spins which couple (predominantly) to a given 
are ideal candidates for such probes, as we shall discuss 
next. 

Nuclear susceptibility. - NMR techniques have been 
widely used to study molecular magnets |ll||l7j]. In the 
following we show that nuclear spins can also be used as 
a local probe to detect coherent tunneling of n. For sim- 
plicity, we restrict our considerations to interactions of 
the form H' = Asi ■ I, which includes both the hyperfine 
contact interaction and the direction-independent part of 
the magnetic dipolar interaction. For 57 Fe, the dominant 
coupling to the electron spin is by a hyperfine contact 
interaction, Ap c s ~ 3.3mK j^J. In contrast, the interac- 
tion of a 1 H nuclear spin with the electron spins is dipolar, 
with a direction-independent term AjSi-I. For AF or- 
der, the sum yields an effective coupling Ajjsn • I, where, 
for Feio, the coefficient Ah — Yli(~ l) i+1 ^4i depends 
strongly on the site of the proton spin I. For many of the 
inequivalent sites, however, Ah is of order 0.1 mK [ p~T|| . 
With N Fc and N R the numbers of NMR-active 57 Fe and 
proton nuclei, as long as -/Vp^Fe + NhAh <C A, the ef- 
fect of the nuclear spins on the electron spin dynamics 
remains small. We define the decoherence rates Ti and 
Ts of the nuclear and electronic spins as the decay rates 
of (Iy(t)I y ) and (n z (t)n z ), respectively. For time scales 
t < 1/Ts, the electron spin produces a coherently oscil- 
lating effective magnetic field with frequency A/H at the 
site of the nucleus. 

In order to show that this field affects the nuclear 
spin dynamics, we now consider a single, NMR-active 
57 Fe or 1 H nucleus (inset of Fig. 1). For k B T < Hujq, 
we may restrict our considerations to the Hilbert space 
spanned by {|<?), |e)}. By using the decomposition = 



with f x = 1 - <D(ir 2 NA/8J) and f y ~ f z ~ 1. It is 
clear from eq. (||) that the susceptibilities Xaa{u) for the 
ring (|l|) do not exhibit resonances at uj = ±A/h. The 
main conclusion we draw from eq. (||) is that a uniform 

magnetic field cannot drive transitions from \g) to |e). (-l) 4+1 sn + S/N of a single spin into staggered magne- 
Starting from the rigid rotor Hamiltonian of the ferric 
wheel 0], this result can also be shown to hold for an 
arbitrary direction of B. 

We consider next the correlation function of a single 



spin, (T T s^ a (T)s ha ) ~ s 2 (T T n a (T)n a ) , for h x < 4Js, 
with i = 1,...,N. In contrast to the correlations of 
S discussed above, (T T Si ia (T)si ia ) indeed exhibits sig- 
natures of coherent tunneling of n. To evaluate the 
correlation function, we use an effective two-state de- 
scription for the ferric wheel and introduce a pseudospin 
Hcr/2, with | f) and | |) being eigenstates of a z . The 
tunneling dynamics of the Neel vector n is then gen- 
erated by the pseudo-Hamiltonian —Aa x /2. Because 
(T T n z (T)n z ) ~ (T T o~ z (r)a z ) in the low-energy sector, we 
obtain immediately 



tization ±sn, n 2 = 1, and fluctuations Sin around the 
Neel ordered state, we obtain H' = As\ ■ I ~ Asm ■ I. 
We show now that, due to (e\n z \g) ^ , the tunnel- 
ing dynamics of n can be obtained from the nuclear spin 
correlation functions (I a (t)I a ). With (e\n z \g) — 0(1), 
the dominant term in H' is Asn z I z , or in the pseudospin 
notation introduced above, 



H' ~ AsnJ z ~ Asa z I z 



(4) 



NMR experiments measure via power absorption the 
imaginary part of the nuclear spin susceptibility, 
X/Qa( w )i an d by pulsed techniques the nuclear spin 
correlation functions (I a (t)I a ) in the time domain p9[ . 
From expanding (I a (t)I a ) in H', 
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cosh[(/3-2|T-r'|)A/2] 
cosh[/3A/2] 
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After analytic continuation, the real-time correlation 
function exhibits the e ±lAt / h time dependence character- 
istic of coherent tunneling. We conclude that local spin 



(I a (t)I a )~(I a (t)I a ) -A 2 s 2 /h 2 dt' At" (5) 

J — OG J — OO 

x ([/,(*'), I a {t)][I z {t"), I a ]) (n z (t>)n z (t"))o, 
it is evident that the dynamics of n enters \" aa ( w )- 
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To evaluate eq. (@), we diagonalize the Hamiltonian 



A 

H = (T T 

2 



liB x I x + Asl z a z 



(6) 



with 77 the nuclear gyromagnetic ratio, which describes 
the ferric wheel in the low-energy sector with a single 
nuclear spin I coupled to Si, eq. (||). We assume ther- 
mal equilibrium for both electron and nuclear spins. For 
Feio m the high field regime, the results may be derived 
by expansion to leading order in As/ A and "/iB x / A, be- 
cause A 3> "fiB Xl As, for both 57 Fe and 1 H nuclei, and 
B x ~ 10T. For I = 1/2, Xi, zz M displays the unper- 
turbed emission and absorption peaks at u> ~ ±.^iB x /h, 
although these are slightly shifted if the hyperfine term 
AS ■ I/N is taken into account. XjxxC^)' however, dis- 
plays resonant absorption and emission of small intensity 
at u = ±(A± jiB x )/H. Finally, 



tanh^^^) fiu- - -.,£,//') (7) 



^) tanh^^(c-A/r,)]-[a 



exhibits satellite resonances at the tunnel splitting A of 
the electron spin system. Their physical origin is readily 
understood in terms of a classical vector model. For A = 
0, l(t) precesses around the static magnetic field B = 
B x e x . For A ^ 0, the coherent tunneling of n leads to an 
additional ac hyperfine field A&i(t) ~ As cos(At/H)e z at 
the site of the nucleus. In contrast to a static hyperfine 
field which induces a change in precession frequency and 
axis, the rapidly oscillating hyperfine field in Feio leads 
only to a small deviation SI(t) = 0(As/A) from the orig- 
inal precession. In particular, 8 I y (t) oc (As/ A) sin(At/fi) 
also oscillates at frequency A/ft and hence gives rise to 
the second term in eq. (Q). 

We restricted the above analysis to the low-energy sec- 
tor of the ferric wheel. To check this approximation, we 
have performed exact numerical diagonalization (ED) on 
a small AF ring with one nuclear spin of / = 1/2 coupled 
to one of the electron spins. For the small systems acces- 
sible by ED, in this case N = 4, s = 3/2, and k z = 0.2 J, 
the field range for which the theoretical framework is ap- 
plicable becomes rather small: 2 J <§; h x <C 6 J. However, 
the numerical results (Fig. 1) for (e|/ y |<?) indicate that 
our analytical value |(e|/j,|<7)| = As/2 A entering eq. (Q) 
is a good approximation. For our parameters, the analyt- 
ical value tends to overestimate (e|/j,|<?) by ~ 30% which 
results from the neglect of n y in eq. (|j). 

We turn next to a discussion of the experimental feasi- 
bility of measuring A by NMR. Because As/ A <C 1, the 
intensity of the satellite peaks at u> = ±A/h (fjj) is small 
compared to that of the main peaks at to = ±"/jB x /h. 
However, this satellite peak intensity may be increased 
significantly by tuning B x close to one of the critical val- 
ues B% at which the magnetization of the ferric wheel 



jumps and A(B X = B^) = 0. Note, however, that our 
theory only applies to high magnetic fields, B x > 7.7T 
for Feio ■ Coherent tunneling of the Neel vector n requires 
A > Ts- From the estimates for the decoherence rate Ts 
given in the introduction we conclude that this condition 
can be satisfied even for a large range of A/Ao -€i 1. 

We consider first 57 Fe, with 7/ = 0A8fi N @. 
For T ~ 2K and B x ~ 10T, the relative inten- 
sity of the satellite peak at A = Ao = 2.18K, is 
(As/A) 2 tanh(/3A/2)/tanh(/3 7/ B a; /2) ~ 0.007. This in- 
tensity, however, increases by a factor of 10 (100) for 
A = 0.1A (A = 0.01A ). For 1 Hwith 7/ = 5.59/zjv @, 
and a typical value As ~ O.lmK, the relative peak inten- 
sity is 2.05 x 10~ 7 (A = A ), 2.25 x 10~ 6 (A = 0.1A ), 
and 2.25 x 10~ 5 (A = 0.01A ). However, the number 
of protons in the ring is much larger than that of NMR- 
active 57 Fe nuclei, 10 < Nn/Np c < 100, depending on 
the doping with 57 Fe. Taking into account that the sen- 
sitivity of proton NMR is larger than that of Fe NMR 
by a factor of 3 x 10 4 |Tj], 57 Fe and proton NMR ap- 
pear to be similarly appropriate for detecting the co- 
herent tunneling of n. The observation of the satellite 
peak in (Q) is feasible, but still remains a challenging 
experimental task. The experiment must be conducted 
with single crystals of Feio (or a Fe6 system with suf- 
ficiently large k z > 2J/(Ns) 2 ) at high, tunable fields 
(10T) and low temperatures (2K). Moreover, because B% 
depends sensitively on the relative orientation of B and 
the easy axis 0,0, careful field sweeps are necessary 
to ensure that A/Ao <gC 1 is maintained | pl| . Note that 
the NMR experiment suggested here could be more easily 
realized with nuclear spins exhibiting higher NMR sensi- 
tivity than 57 Fe. 

We show now that, from NMR spectra, also an upper 
bound for Ts can be extracted. The NMR resonance 
lines are broadened by the decoherence of the nuclear 
spin, with width Tj at ±jjB x /h. The NMR resonances 
at oj = ±A/h also involve correlation functions of n, see 
(||). Thus the decoherence of the electron spin, Ts, adds 
to the linewidth, and the width of the satellite peak, 8, 
is bounded by T/ + < S. Measurement of 8 then pro- 
vides an upper bound for Ts- Further, Ts — A marks the 
transition from coherent to incoherent tunneling dynam- 
ics. Hence, if S < A this would indicate unambigously 
that quantum tunneling of n is coherent. Note that the 
maximum peak height of a Lorentzian resonance line is 
0(1/(5), so a large Ts (< A) would make detection of the 
satellite peak increasingly difficult. 

ESR measurements in the presence oj hyperfine in- 
teraction. - We show now that, in the presence 
of a hyperfine coupling H' , the electron spin sus- 
ceptibility of the ring, Xaa(u), also exhibits the sig- 
nature of a coherent tunneling of n. This results 
from the fact that integration over the initial and fi- 
nal nuclear spin configurations causes the matrix ele- 
ments (e|S|g) occurring in the spectral representation 
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of Xa«( w ) ( ec l- (@)) to become finite. In the high 
field limit h x 3> tvjJo, the matrix elements become 
Kel^lp)! ~ (As/2h x )(A S/Wk z s 2 )\cos{irNh x /M)\ 
and |(e|5 z |g)| ~ As/2h x . For a small system (N = 4, 
s = 3/2), we have again confirmed the qualitative fea- 
tures of these results by ED. It follows that Xaa( w ~ 
A/h) = 7r|(e|S Q |5)| 2 tanh(/3A/2)(5(w- A/h) exhibits res- 
onances at A/h. A qualitative understanding of this 
result may be obtained by noting that a nuclear spin 
polarized along B results in an effective magnetic field 
e x A/2gfj,B acting only on si. In a classical descrip- 
tion this hyperfine field causes S to acquire a component 
nAs/2h x along n, and the coherent tunneling of n now 
also results in an oscillation of the total spin S. Again, 
the decoherence rate of these oscillations, and hence the 
linewidth of the ESR resonance, is bounded by T$- Note 
that As/2h x = 1.2 x 1(T 4 for a single 57 Fe nucleus with 
B x = 10T, so the ESR signal is very weak in this case. 

However, our calculations apply to any impurity spin 
j, which interacts with a single electron spin only, H' = 
As\ • j. In particular, for an electronic j, A is typically 
10 3 times larger than for a nuclear spin, and ESR tech- 
niques become a valuable tool for detecting the tunnel- 
ing of n. One advantage of this technique is that it 
is no longer necessary to have A/Aq < 1 to obtain a 
large signal intensity, and thus the complete range of 
tunnel frequencies could be explored experimentally . 
Our calculations also apply to situations in which sev- 
eral impurity spins produce different net magnetic 
fields for the two sublattices of the AF ring. For illus- 
tration, we discuss two simple scenarios. We consider 
N/2 impurity spins ji (jj = 1/2) coupled to si, s 3 , . . . , 

H ' = A 12ili J2i-1-S2i-1- For h x > As, k B T, all j, align 
with the magnetic field B. Since they all couple to one 
sublattice only, their net magnetic fields acting on add 
up, \(e\S z \g)\ ~ (N/2)As/2h x , leading to a (iV/2) 2 -fold 
enhancement of the ESR-signal of a single impurity. In 
contrast, a single impurity j coupled to both si and S2, 
H' = A] ■ (si + S2), results in the same net magnetic field 
acting on both sublattices, (e\S z \g) = 0. 

Conclusions - We have shown that NMR and ESR 
techniques can be used to measure both the tunnel split- 
ting A and the decoherence rate T$ in the ferric wheel. 
For Feio, we showed that our proposal is within experi- 
mental reach. Our considerations apply to any AF ring 
system described by Ho, eq. ([!]), with some impurity 
spin coupled to one of the electron spins. Hence, the 
proposed scheme may prove useful for a wide class of 
molecular magnets. 
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FIG. 1. Comparison of the analytical result for A (upper 
panel) and (e\I y \g) (lower panel) with exact numerical re- 
sults for a small system, N = 4, s = 3/2, k z — 0.2J, and 
A — 9 x 1CP 5 J. The numerical values (o) for (e\I y \g) are well 
approximated by As/2 A (solid line), with a small offset in 
B x resulting from the shift of the magnetization steps when 
k z 7^ 0. For reference, the ratios As/2A with the numerical 
values for A are shown (•) in the lower panel . 



